12. Sept. 1991 



Isospin Dependent Effective Interaction 



in 

On 



Q 

oo 



in 

OS 



o 

=5 



X 



in 



Nucleon-Nucleus Scattering 



Taksu Cheon 
Department of Physics, Hosei University 



Chiyoda-ku Fujimi, Tokyo 102, Japan 



and 



> 

en 

q . Kazuo Takayanagi 



Department of Mathematical Sciences, Tokyo Denki University 
Hatoyama, Saitama 350-03, Japan 



Abstract 

The isospin-dependent component of the effective nucleon-nucleon interaction that 
causes the AT = 1 (p, p') and (p, n) reactions off nuclei is studied. It is shown 
that the corrections to the impulse approximation comes from the (7-matrix type 
correction and the rearrangement term. They are numerically estimated with the 
isospin-asymmetric nuclear-matter reaction matrix approach. The analysis of the 
isobaric analog transitions 42 Ca(p,n) 42 Sc and 48 Ca(p,n) 48 Sc are presented. 
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The nucleon direct reaction of few hundred MeV to GeV projectile energy 
has become one of the preferred tools to extract the nuclear structure informa- 
tion complementing the reaction with the electromagnetic probes. Its success 
crucially depends on the good command of the effective nucleon-nucleon inter- 
action used in the calculations based on the distorted-wave Born approximation 
(DWBA). From the theoretical point of view, there is a possibility of performing 
a convergent calculation of Watson expansion [1,2] or its variants [3], since the 
impulse-approximation is expected to be a reliable starting point at these ener- 
gies. In practical analysis with the currently available computational resources, 
however, the treatment based on the density-dependent effective interaction [4,5] 
seems to be the only workable scheme to handle the medium corrections. 

For the elastic scattering, the optical model potential obtained by folding 
the Brueckner reaction matrix (g-matrix) [6] is known to represent the leading 
terms in the Watson expansion [7], and the g- matrix approach has been reason- 
ably successful in describing the elastic scattering experiments [4,5,8,9]. On the 
contrary, for the inelastic scattering, there had been no such general recipe. Many 
DWBA calculations had been performed, however, with the density dependent 
(/-matrix without any theoretical justification. In this situation, we have shown 
for the simplest case, i.e. the isoscalar natural parity transitions, that the tran- 
sition potential is the sum of the ^-matrix and the rearrangement term which 
is written in terms of the derivative of the ^-matrix with respect to the nuclear 
matter density [10,11]. Attempts to describe experimental data along this line 
have been pursued in recent years, and met a great success. The set of density- 
dependent effective interactions constructed by Kelly and collaborators [12-14] has 
been proven quantitatively accurate and practically useful to describe the elastic 
scattering and isoscalar natural parity excitations. In this background, we believe 
it quite urgent to identify medium corrections in other components of the effective 
nucleon-nucleon interaction than the spin-isospin independent part, since one of 
the advantage of the nucleon projectile over the electron is the wider variety of 
nuclear excitation modes accessible by it. 

Although the isobaric analog excitation [15] has been one of the best-studied 
nuclear excitation modes, the isospin-dependent part of the effective interaction 
has been given very little theoretical attention. This unsatisfactory state of affairs 
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needs to be altered before any systematic empirical analysis be performed. In 
this Letter, we intend to identify the leading medium corrections to the isospin 
dependent component of the effective nucleon-nucleon interaction which induces 
the isovector nuclear excitations. In the formulation of such a process, the g-matrix 
in an asymmetric nuclear matter where the Fermi momentum of the proton differs 
from that of the neutron naturally appears. This (7-matrix is used to evaluate the 
medium effects on the isovector transitions. We report the result of a numerical 
example at the incident energy E p = 150MeV. It is applied to the isobaric analog 
transitions 42 Ca(p,n) 42 Sc and 48 Ca(p,n) 48 Sc at E p = 135MeV. We hope that this 
work is to be the first step for the microscopic analysis of all the components of 
effective nucleon-nucleon interaction for intermediate energy nucleon scattering. 

The (7-matrix in an isospin-asymmetric nuclear matter [16] satisfies the Bethe- 
Goldstone equation formally identical to the one in the symmetric nuclear matter: 

g = v + vGoQg , (1) 

where v is the free nucleon-nucleon interaction, and Go and Q are the two-nucleon 
propagator and the Pauli operator in the medium, respectively. We neglect the 
self-energy term in Go assuming that the incident proton has sufficiently high en- 
ergy. Then the difference of proton and neutron densities, p p and p n , generates the 
isovector component of Q, in addition to the isoscalar component which charac- 
terizes the isospin-symmetric nuclear matter. As a result, the effective interaction 
g of eq.(l) has four independent components g T ,r z specified by the total isospin T 
and its ^-component T z of the two interacting nucleons. We write the (7-matrix in 
the following way, featuring the tensor properties in isospin space: 

g = g [o] +g M Tl . T2 + g H {Tl + T2)+g [P] y ^ [TlXT2] (2) ^ (2) 

where r\ and T2 are the isospin operators of two interacting nucleons. The four 
coefficients in the eq.(2) are written in terms of g T ,r z as 



' [0] = ^ (dii +9io + 9i-i +9oo) , 



(3a) 



9 [t] = ( gi1 +9l ° + 9i-i ~ 3#oo) , (36) 

g [a] = \ (911-91-1) , (3c) 

9 [P] = \(9n- 2^io + <7i-i) • (3d) 



From eqs.(3a) to (3d), one recognizes that g^°\ and are invariant with 
respect to the exchange of protons and neutrons in the medium, i.e. p p <-> p n , 
while g^ changes its sign. These properties put restrictions of the possible form 
of density dependence of each term. A transparent view is obtained by defining 
the isoscalar density p 1 ^ = p p + p n , and the isovector density p^ = p v — p n . 
Specifically at the isospin-symmetric limit p^ =0, one obtains 

Also, taking the leading density-dependence of g^ a \ one obtains 

We now consider the elastic and inelastic scatterings in order. The elastic 
optical model potential U of a nucleon in the isospin asymmetric nuclear matter 
is given by folding the effective interaction g with the asymmetric nuclear density. 
It is convenient to split the resultant optical model potential into isoscalar and 
isovector part d la Lane as [17] 

U = U [s] + U [v] t z , (5) 

where each term is given by 

uW=gMpM+gWpM 



g [0] p [s] 



(6a) 



and 
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U M =g [ a ] p [s] + ^ g [r] +g [^ p [v] 

-( PW ^ [Q1+ ^) PM • 

In the derivation of t/"^ and U^ v \ we have dropped all terms that are of second 
or higher orders in powers of p^ v \ In the last line of eq.(6b), the relation (4b) is 
used. The above expression shows that the optical model potential U of eq.(5) is 
obtained by folding 

Va =gl°l p M + frW° gM+gM) Tl .^ 

dpM ( 7 ) 

= V°l + V ll r l ' r 2 > 

with the nuclear densities p^ and p^. One can regard v e \ defined by eq.(7) as the 
effective interaction for the elastic scattering. 

For the inelastic scatterings, both through the diagrammatical analysis of 
Watson expansion and the macroscopic collective excitation model, one can relate 
the inelastic transition potential and the density derivative of the optical potential 
U for the elastic scattering [10,11]. Assuming that the transition of the target 
nucleus is described by the isoscalar and isovector transition densities p\"J and py , 
one can write the transition potential Ut r which describe the inelastic scattering 
as 

This yields, for N = Z nuclei (pM = 0), the expression 

which explictly shows that the transition potential can be obtained from the ef- 
fective interaction 



Vtn = V l S+V^ Tl -T2 , (10) 
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where the isospin-independent and dependent components, ujj and v^, are given 
by 

^ 1= * |0|+ " W ^ 101 ■ f 11 "' 

and 

The second term of eq.(lla) is the so-called rearrangement term for the isoscalar 
transition [10] which roughly doubles the density dependence of the inelastic ef- 
fective interaction. Eq.(llb) is one of the main result of this Letter. Its second 
term is analogous to the one in eq.(lla), and no less remarkable since it persists 
at p[ v i = despite that itself vanishes at this limit. We now compare eq.(7) 
with eq.(ll), i.e. the effective interactions for elastic and inelastic scatterings. For 
the isospin-independent component, one recognizes 

«S = ( 1+ " W ^)^ ■ < 12 > 

which has been instrumental in determining the empirical density-dependent in- 
teraction [12-14]. For the isospin-dependent component, one finds 

V M _ Jr\ _ Jt] , o [8]_0_ a l<*] ( 13) 
V in~ V el -9 +P Q p[v] 9 ■ K±0) 

This shows that the extra factor appearing in eq.(12) can be formally eliminated 
in the case of isospin-dependent component. We stress that this does not mean the 
absence of the rearrangement term. It is achieved by natural but purely formal 
redefinition, eq.(7). In principle, therefore, with eq.(13), the isospin-dependent 
empirical interactions can be constructed from the combined analysis of the elastic 
scattering and the AT = 1 inelastic and charge-exchange reactions on N ~ Z 
nuclei. 

We now turn to the numerical assessment of the medium correction to the 
isospin-dependent component of the effective interaction. The detailed description 
of the technique to solve eq.(l) in an isospin-asymmetric nuclear matter and the 
full results will be published elsewhere. Here, we show the relevant result of an 
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example of such calculation. At the projectile energy of E p = 150MeV, Reid Soft 
Core potential as the input free two nucleon potential, we obtain the following 
number for the volume integrals of isospin dependent components of ^-matrix. 

g [r] = [{22 + 75i} + {-24 - 53i}p [s] ] + [{8 + 22i} + {17 - 17i}p [s] ]L • S , (14) 
g [«l = [{6 + 380p H ] + [{-1 + 10p H ]L-S , (15) 

where the densities p 1 ^ and p^ are normalized to the nuclear matter saturation 
density. The isospin dependent component of the ^-matrix, g^ T \ is identical to the 
one in the symmetric nuclear matter in ref. [18]. It includes the density dependence 
induced by the nonlocality. Only the contribution from the allowed channel is 
included, and thereby the antisymmetrization is taken into account. The numbers 
in eqs.(14) and (15) can be also interpreted as the strength of the interaction 
in zero-range approxiamation. One observes that, at this energy, the density 
dependence of tends to cancel that of g^ T \ which already exists in the g- 
matrix in the symmetric nuclear matter, in the calculation of the isospin-dependent 
effective interaction, eq.(13), especially in the imaginary part of the central force. 
More to the point, these two density-dependent corrections are of the same order 
numerically. It clearly shows that the partial inclusion of the medium correction, 
only with (7-matrix, for example, cannot be very meaningful. 

In order to show how the density-dependent effective interaction eqs. (13-14) 
works in the inelastic scattering, we show the result of DWBA calculations [19] for 
the isobaric analog transitions 42 Ca(p,n) 42 Sc and 48 Ca(p,n) 48 Sc at E p = 135MeV. 
For the facility of the calculation, we adopt following two approximations: First, 
we take the empirical optical model potential by Schwandt et al [20], instead of 
calculating it microscopically. Second, we replace the density-independent part of 
eqs. (14) and (15) by the three-range Yukawa parametrization of free two-nucleon 
scattering matrix by Franey and Love [21] at E p = lAQMeV (FL140). The density- 
dependent portion of the above effective interactions is treated as Yukawa force 
of very short range. The calculated differential cross section is shown in Figures, 
where dashed line corresponds to the impulse approximation with FL140 interac- 
tion; dotted-dashed line, the result with of eq.(llb) being replaced by 
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solid line, the result with the full medium correction included. We can see the 
following effect: First, the calculation with the t-matrix overestimates the experi- 
mental cross section [22]. Second, by changing the interaction to the g-matrix, the 
calculation lowers the cross section, but overshoot the data. Last, by including all 
the medium effects, the calculation hits the experimental data points very well. 
A warning is due, however, to the overemphasis on the favourable comparison 
with the experiment at this point, since the uncertainty of the isospin-dependent 
component of the free t-matrix is already very large. Also, the approximations 
in the current calculation leaves some space for further improvement. Rather, we 
would reiterate our basic point that the medium modification should be, and can 
be estimated with the proper inclusion of the rearrangement term. 

We acknowledge our gratitude to Prof. K. Yazaki for the valuable discussions. 
Thanks are also due to Prof. H. Sakai for his advice in experimental and compu- 
tational aspect of this work. The numerical work has been performed on SPARC- 
Station2 at Hosei University, Department of Physics, and on VAX 6000/440 of 
Meson Science Laboratory of University of Tokyo. This work is financially sup- 
ported by Research Institute for Technology, Tokyo Denki University. 
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Figure Captions 



Fig. [I]: The differential cross section for the reaction 42 Ca(p,n) 42 Sc. The dashed 
line is the result of the impulse approximation. The dotted-dashed curve partially 
accounts for the medium corrections, while the solid line includes all the medium 
modifications. 

Fig. [II]: The differential cross section for the reaction 48 Ca(p,n) 48 Sc. Notation is 
the same as for fig. [I]. 
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Fig. 1 




Fig. 2 



